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The canonical quantization of a Schwarzschild black hole yields a picture of the 
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tropy of the black hole by explicitly counting the microstates corresponding to a 
macrostate of fixed total energy. 
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1. Introduction 

It has been recognized for some time that black holes behave as thermodynamic 
objects with a characteristic temperature and entropy and that these quantities 
are inherently quantum mechanical in nature. 1 ' 2 This makes a clear understanding 
of the origins of black hole thermodynamics in terms of statistical principles one of 
the more interesting open problems, because a microscopic description of the black 
hole entropy requires a quantum theory of gravity so that detailed first-principle 
investigations concerning the black hole entropy should contribute toward a better 
understanding of how such a theory may be constructed and interpreted. The 
earliest attempt at a microscopic theory of black holes was due to Bekenstein 3 , 
who concluded that the horizon area is the analog of an adiabatic invariant in 
mechanics. He then invoked the Christodoulou-Ruffini process 4 and the Bohr- 
Sommerfeld quantization rules to argue that the eigenvalues of the horizon area 
operator, A, of the black hole must be equally spaced, A n ~ nip. Dividing the 
horizon in cells of Planck area which get added one at a time, and assuming that 
each cell has the same (small) number of states, say k, he was able to derive the 
area law of black hole thermodynamics by estimating the number of microscopic 
states to be f2 « k n . 

Bekenstein's hypothesis has led to many interesting attempts to derive the area 
quantization law and the black hole entropy from first principles. These attempts 
have mostly taken either the loop (or canonical) quantum gravity 5 approach or 
the string theory 6 approach. The ability to reproduce the Bekenstein-Hawking 
entropy from a microscopic counting of states must be considered a measure of the 
success of a candidate quantum theory of gravity. Thus, the successful derivation of 
the Bekenstein-Hawking entropy, for example, from a microcanonical ensemble of 
D— brane states 7 is generally considered a triumph for string theory and a necessary 
condition for string theory to be a convincing candidate for a theory of quantum 
gravity. 

An approach to the entropy problem that makes explicit reference neither to 
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string theory nor to canonical quantum gravity is to count the number of states 
for a conformal field theory corresponding to the asymptotic symmetries of the 
black hole or to the symmetries of its horizon. This mechanism for calculating the 
number of microstates was recently proposed by Strominger 8 for the BTZ black 
hole. 9 Strominger used a result of Brown and Henneaux 10 which says that the 
asymptotic symmetry group of AdS3 is generated by two copies of the Virasoro 
algebra with central charge 3Z/2G, and combined this result with Cardy's formula 11 
to compute the asymptotic growth of states for this conformal field theory. A 
semi-classical analysis of a hot black hole suggests that the horizon is thermally 
oscillating while maintaining a fixed area. Recently Carlip 12 has argued that 
Strominger's result is more generic than was originally believed, because the algebra 
of surface deformations of any black hole in any dimension contains a Virasoro 
algebra consisting of deformations that leave the horizon fixed. Cardy's formula 
once again yields the correct Bekenstein-Hawking entropy. In this approach, the 
horizon is treated as a boundary and all of the relevant degrees of freedom of 
the black hole are assumed to lie on it. The states themselves are not explicitly 
displayed. 

The results from loop quantum gravity have been no less dramatic. Here spin 
network states are used as a complete, orthonormal basis in the Hilbert space and 
loop microstates are explicitly counted using techniques developed in ref.[13]. In 
one approach 14 , loop quantum gravity is used to compute the microstates of the 
horizon, holding the area fixed. In another approach, one analyzes the classical the- 
ory outside, treating the horizon as a boundary. Quantization yields surface states 
which are counted by an effective Chern-Simons theory on the boundary. 15 Both 
approaches yield an entropy proportional to the area of the horizon, but the pro- 
portionality constant is a free, finite and dimensionless parameter, not determined 
by the theory and the horizon area eigenvalues are not equally spaced. 

In this article we take a midi-superspace approach to the canonical quantization 
of the Schwarzschild black hole and show that it leads naturally to Bekenstein's 
mass quantization law (equally spaced area eigenvalues) as well as to the area law 
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of black hole entropy as computed from a genuine microcanonical ensemble. The 
quantization leads to an amusing picture of a black hole which is akin to statistical 
bootstrap models for hadrons, whose statistical mechanics was studied many years 
ago by Frautschi 16 and Carlitz 17 . We will use similar techniques to study the 
statistical properties of the Schwarzschild black hole, but our starting point will 
be a recent solution 18 of the Wheeler-DeWitt equation in terms of variables first 
introduced by Kuchar and Brown. 19 ' 20 



2. Quantization 

As we have previously shown 18 , combining the Hamiltonian reduction of spher- 
ical geometries due to Kuchar 19 and the coupling to dust as proposed by Brown 
and Kuchar 20 allows for the derivation of a simple, decoupled (Wheeler-DeWitt) 
equation describing the Schwarzschild black hole. The non-rotating dust is intro- 
duced in such a way that its role is only as a time keeper. Consider the gravity-dust 
system, 

S = --1-J d 4 x^TZ - l -J d 4 x^e(x) [g a(3 U a uP + l] , (2.1) 

in the general spherically symmetric spacetime, 

ds 2 = N 2 dt 2 - L 2 {dr-N r dt) 2 - R 2 dQ 2 , (2.2) 

where N(t,r) and N r (t,r) are respectively the lapse and shift functions, R(t,r) is 
the physical radius or curvature coordinate, e(t,r) is the density of the collapsing 
dust in its proper frame, 1Z is the scalar curvature and U a are the components of 
the dust velocity. The action, S, in (2.1), may be be recast into the form 



dtdr 



P L L + P R R + P T f - NH - N r H r 
+ surface terms, 



(2.3) 



where we have introduced, following Brown and Kuchar, 20 the dust proper time 
variable, r, which in general will serve as extrinsic time. Pl and Pr are the 
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momenta conjugate to L and R respectively, and the super Hamiltonian, H, and 
super momentum, H r , are respectively given by 



H 



PlPr 
R 



LPl 
2R 2 



+ 



L 
2 



FT 
~2L 



+ 



RR' 



+ P T yfl + r' 2 /L 2 



(2.4) 



and 



H r = R'Pr - LP' L + r'P T , 



(2.5) 



where the prime denotes a derivative with respect to the ADM label coordinate 
r. The constraints in the above form do not "decouple" and are very difficult to 
resolve as they stand. However, from the general system in (2.2), Kuchaf 19 showed 
how one can pass by a canonical transformation to a new canonical chart with 
coordinates M and R together with their conjugate momenta, Pm and Pr, where 
M is the Schwarzschild "mass" and R is the curvature coordinate. In this system 
the constraints are greatly simplified and the phase space variables have immediate 
physical significance. The canonical transformation is well-defined as long as the 
metric obeys standard fall-off conditions 19 and, as long as these fall-off conditions 
are obeyed, the surface action can be recast in the form 



surface terms = J dt [7r + f + + 7r_f_ — N + C + — iV_C_] , 



(2.6) 



where r± are the proper times measured on the parametrization clocks at right 
(left) infinity. The constraints C± = ±7r± + M± identify their conjugate momenta 
as the mass at right (left) infinity. In terms of the new variables the entire action, 
along with the surface term is 



S = 



dt / dr 



P M M + P R R + P T f - NH - N r H r 



(2.7) 



where Pm = Pm — t' , P t = P t + M' and Pm and P T are the original Brown- 
Kuchar variables. The transformations leading up to these variables may be found 
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in ref. [19]. In passing to the transformed momenta, P T and Pm-, we have made a 
canonical transformation 18 generated by Mr', which effectively absorbs the surface 
terms. (Thus, we have implicitly fixed the dust proper time to coincide at infinity 
with the parametrization clocks.) 

The super Hamiltonian and super momentum constraints become 



H = 



F- l M'R' + FP R (P M + r 1 ) 



+ {P T -M')^l + r' 2 /L 2 = 



(2.8) 



and 



H r = M'P M + R'Pr + r'P T = 0, (2.9) 
where we have used 



L = F-'R' 1 - F(P M + r'Y (2.10) 

and F — 1 — 2M/R. L 2 , being the component g rr of the spherically symmetric 
metric in (2.2), must be positive definite everywhere. F is positive in the exterior 
(Schwarzschild) region and negative in the interior and this will play an important 
role in the consistency conditions that follow. By direct computation of Poisson 
brackets, it is easy to determine the "velocities" in terms of the conjugate momenta 
from the above expressions and they are 



f = N^l + r' 2 /L 2 + N r r', 

a = -wcw) + (2 . n) 

Lj 

NR'r'(P T -M') . 

M = ^4 + N r M'. 

L 6 

The constraints in (2.8) and (2.9) generate the local symmetries of the theory. In 
the Schroedinger representation, they correspond to functional differential equa- 
tions that the state functional must satisfy, i.e., the canonical variables must be 
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raised to operator status and the constraints considered as operator constraints on 
the (Wheeler-DeWitt) wave functional i/?wd[t, R, M], 

Hij WD [r,R,M] = = H r il> W D[T,R,M}. (2.12) 

These equations are, in fact, an infinite set of equations, one for each spatial point 
on the spatial hypersurface. 

If the super momentum constraint in (2.9) is used to eliminate Pm m the 
expression for the super Hamiltonian in (2.8), the latter constraint turns into 18 

(Pr-M'f + FP 2 R - = 0. (2.13) 

We will now specialize to the black hole by requiring M' = so that only the 
homogeneous mode of M(t, r) survives. In this way the dust is made tenuous and 
the proper time variable must be thought of as the proper time of a test particle 
(an ideal clock) in free fall. To understand its role in the canonical reduction, recall 
that the Schwarzschild metric in comoving coordinates can be viewed as a special 
case of a marginally bound Tolman-Bondi metric for the collapse of inhomogeneous 
dust, which has the general form 

ds 2 = dr 2 - R' 2 (T,p)dp 2 - R 2 {T,p)dQ 2 , (2.14) 

where r is the dust proper time and the prime denotes a derivative with respect to 
p. The curvature coordinate, R(r,p), is obtained from Einstein's equations, after 
an appropriate scaling 21 , as 

R 3/2 = P 3/2 " \^Fjp)T, (2.15) 

in terms of an arbitrary function (the "mass" function), F(p). It represents the 
mass contained within a shell of radius p. The energy density of the dust is given 
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by 



= (2.16) 



In particular, if -F(p) = a 2 (const.), the energy density vanishes everywhere except 
at the singularity, R = 0, and the metric in (2.14) describes a (Schwarzschild) black 
hole of mass a 2 /2. This follows from the coordinate transformation 22 



T + a 



,3/2 = R 3/2 + + a J 



VRdR 

(2.17) 

RdR x 



R-a 



2 I ' 



which takes (2.14) into the standard Schwarzschild form, with 2M = a 2 and where 
T is the Killing time. This spacetime transformation may be re-expressed as a 
point transformation on the phase space. 

The resulting Wheeler-DeWitt equation (with M' = 0) is decoupled, and the 
quantum state may be expressed formally as a direct product of states defined at 
each spatial point, 

w>wd) = n ( 2 - ig ) 

r 

where we have used r as a label. Each of the \^ r ) is normalized w.r.t. a suitable 
measure in the space of functions r(r), R(r), M at fixed r and the system reduces 
to a set of independent Schroedinger equations, one for each spatial point. The 
wave functional is given by 

ij WD [r,R,M} = T[* r [r(r),R{r),M\ 

(2.19) 

= #i ® * 2 ® ••• ® ••• <8> ^N- 

The second equation is written by imagining that a lattice is placed on each spatial 
hypersurface so that the classically continuous label coordinate, r, is discretized. 
In this form, the Wheeler-DeWitt wave functional represents a collection of, say, N 
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decoupled systems, each determined by the same Schroedinger equation and obey- 
ing the same boundary conditions. The precise value of iV cannot be determined 
at this level. We will ascertain its value by requiring it to maximize the density of 
states. 

The second constraint enforces spatial diffeomorphism invariance of the wave 
functional on hypersurfaces orthogonal to the dust proper time. By taking func- 
tional derivatives of ipwD m (2-19), it is easy to see that this constraint enforces 
\&J.[t(t"), R(r), M] = 0, where the prime denotes a derivative with respect to the 
label coordinate r. 

We also see, from eq. (2.13), that, for every label r, the (Schroedinger) equation 
reads 

V 2 ^ = 7 a6 V a V^ = HV = 0, (2.20) 

where r y a i > is the field space metric, r y a i ) = diag(l, 1/F), and V a is the covariant 
derivative with respect to this metric. The operator V 2 is the Laplace-Beltrami 
operator having 7^ as its covariant metric, and eq. (2.20) is a massless "Klein- 
Gordon" equation. It is hyperbolic in the region R < 2M (the interior of the 
Kruskal manifold) but elliptic in the region R > 2M (the exterior). This is because 
the quantity F is negative in the interior, but positive in the exterior. As shown 
in ref. [18], this means that the unique positive energy solution of the Wheeler- 
DeWitt equation in the exterior, that is compatible with spatial diffeomorphism 
invariance, is identically zero. The dynamics is therefore confined to the interior of 
the hole. This is consistent with the assumed geometry, as the asymptotic observer 
sees the exterior region of the spacetime as static. 

In the interior, the Schroedinger equation is hyperbolic and it is convenient to 
transform to the coordinate R* defined by 

R* = - ^R{2M - R) + Mtan" 1 



R-M 



^R(2M - R) 



(2.21) 
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The new coordinate lies in the range (—^-, + 1L ljr-) and the wave equation, 

d 2 T ^ - = 0, (2.22) 

now defines the quantum theory whose Hilbert space is Ti := £ 2 (R, dR*) with 
inner product 

I 7rM 

2 

(*i,*2> = y di2„#{tf 2 • (2.23) 



7T-M 

2 



The general (positive energy) solution is 

= c + (M)e- iE ^ + ^ + c_(M)e- li5 ( r ^*) (2.24) 

where c± are functions only of M. We must impose the super momentum con- 
straint, which reads 

(t' + R')c+(M)e- iE ( T+71 ^ + (r / - J R / )c-(M)e-^ r -^) = 0, (2.25) 

assuming E > 0. A consistent and physically meaningful solution to this equation 
is t' = = 0. Returning to (2.11), we see that the choice implies that f = N and 
M = 0. Setting N — 1, the dust proper time turns into the asymptotic Minkowski 
time and the energy, E, should be associated with the ADM mass of the black 
hole. 

Imposing continuity across the horizon, this solution will match the solution 
in the exterior = 0) at R = 2M, if 

c_(M) = - c+(M)e- tEMn , (2.26) 
so that the solution in the interior is now of the form 18 

^> in = C+(M) L-iE{T+H) _ e -iEMn e -iE(r-R^ _ 

There does not seem to be a natural way to impose further boundary conditions. 
Certainly, boundary conditions cannot be imposed at the classical singularity where 
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the canonical reduction will break down anyway. Nevertheless, one notes that the 
parity operator R* — > — R* commutes with the Hamiltonian. States of definite 
parity will vanish at the classical singularity and exhibit a discrete spectrum given 
by 

= -^=e- lET cos ER„ EM = (2n + 1) , 

VnM (22Q) 

y ( r ] = -^=e~ iEr sin ER* EM = In. 

VirM 



" in 



Furthermore, as the dust proper time is identified with the asymptotic Minkowski 
time and the total energy with the ADM mass of the black hole, we are led to the 
Bekenstein mass quantization rule 

M n = ^fnM p (2.30) 

for the definite parity states. 

States of indefinite parity do not admit a quantized mass spectrum, but there 
are three good reasons to confine attention to states of definite parity. Firstly, 
because the parity operator commutes with the hamiltonian, states of definite 
parity are guaranteed to remain so for all time, which is in harmony with our 
intuitive notion of an "eternal" black hole. Secondly, definite parity eigenstates do 
not support the singularity at the origin, which, given that the entire canonical 
quantization program breaks down there, is an attractive feature. Finally, and 
perhaps most importantly, we shall count only the definite parity states in what 
follows and show that they fully account for the entropy of the black hole. 

3. The Entropy 

To compute the entropy, we must enumerate the states of the system. To this 
effect, it is convenient to reformulate the problem by recognizing that the wave 
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equation at each label, r, in (2.20) is derivable from the action, 



1 N r 

S = d 2 X^\r b d a ^ r d b ^ r , (3.1) 



K 



where X e (r, i?*), the integral is over the interior of the Kruskal manifold and 
such that * r (r,!R*) = at = -ttGM/2,+7tGM/2, V r. Imposition of these 
boundary conditions automatically confines attention to states of definite parity. 
Recall that the Wheeler-DeWitt wave functional is a direct product state, as given 
in (2.19). As we will see, each component of the direct product can be thought of 
as describing a tower of oscillators living in the internal space parametrized by the 
phase space coordinates (r, R*). 

Using (3.1), performing a mode expansion of \EV and combining both even and 
odd parities, we can express the contribution of any one of the lattice sites to the 
total energy of the system in terms of pairs a n , f3 n of creation and annihilation 
operators as follows 

M 2 „ , 

H r = ^J2( a l a -r + A r Pn r ), (3.2) 

where 

K r ,o4j = n r 

(3.3) 

[Ai r ,/4j = n r 

are the only non-vanishing commutators. At each site, r, one therefore has a hier- 
archy of two dimensional oscillators. The total energy is the sum over contributions 
from each of the iV sites, i.e., 



M 2 N 

r=l n r 
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which gives (the total energy is the mass of the black hole) 



1 



M 



M 



r 



(3.5) 



y 1 T*bj* 



M 



vWm p , M e Nu{0}, 



where n r ,l r and N Ur ,Ki r are respectively the level number and the occupation 
number at level number n r (l r ), corresponding to the oscillators at site r. 

Let pd(t?7t) be the density of levels describing each site, r. This is just the 
number of states with mass given by m 2 = v r M^ {y r £ N), and is known to have 
the asymptotic (m r >> M p ) form 23 



where c is a constant and D is the dimension of the oscillator. For a generic level 
density, p/)(m r ), the density of states may be written as 



where Mo is the lowest value of m r for which the density of levels is valid. As 
each site is fixed, there are no further phase space integrals. The delta function 
in (3.7) imposes energy conservation as required by (3.5). Let us assume that 
PD( m r) is such that the dominant contribution to the mass integrals comes from 
states with large mass (unless these states are forbidden by energy conservation). 
This is certainly true for the level density in (3.6). Then define the quantity 




(3.6) 




(3.7) 




(3.8) 
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in terms of which we may write 

m,M) = ^n(M). (3.9) 

The S— function restricts the limits of the mass integrals in the definition of cfn(M), 
they no longer run to infinity. Following Carlitz 17 , we estimate the r th integral by 



i 

/ 



v /MA r (a 

dm r p]j(m r ), (3.10) 

Mo 

provided that the A r (x) are subject to the constraint 

N 

J2 A r(%) = X. (3.11) 

r=l 

The maximum contribution to aj^(M) is obtained when the A r (x) are all of the 
order of x/N. This provides an estimate for the integrals in (3.8), and one finds, 
quite generally, 

n(N,M) = a N f(0 N , (3.12) 

where £ = bM 2 /NMp, and a and b are dimensionless constants . The most probable 
number of sites is obtained by maximizing Q(N, M) with respect to N. Therefore, 
consider 

^\nQ(N,M) = In a + In/(0 - ^ln/(0 = 0, (3.13) 

which depends only on ^. Assume that f2 is maximized for some value, say a -1 , 
of £. Then this value gives the corresponding number of sites as 

M 2 

Nmax = db—^. (3.14) 

IVlp 

Nmax is therefore proportional to the area of the horizon. A remarkable conse- 
quence is that all the degrees of freedom can be treated as though they resided just 
there, that is on the horizon itself. 
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Equation (3.13) can be integrated and the solution written in terms of a single 
dimensionless constant, 7. One finds that 

e-*/(0 = -■ (3-15) 

(X 

This constant, 7 then determines the maximum number of states according to 

\nQ max (M) = N\n(a) + N\nf(C)\ Nmax 

= Nlne-^m + AHn(a) + N^\ Nmax ^ 

M 2 

= N^\ Nmax = 67—5 

ivip 

and is, itself, to be determined from (3.15) at £ = a" 1 , i.e., 

7 = aln[a/(a -1 )]. (3.17) 

One thus recovers the entropy 

S = \nn max (M) = 4^2 (3-18) 

where A is the horizon area. This is the area law, provided that 7 > 0, and it is 
independent of the precise form of the density of levels (except for proportionality 
constants), only requiring that the latter is such that the dominant contributions 
come from the most massive states permitted by energy conservation. S inherits 
its dependence on M only from the dispersion relation in (3.5). 

As an example, let us use the density of levels as given in (3.6) although we 
will see that it is strictly not correct to do so. In the present context D = 2, giving 



P2( m r) = c x m r 3//2 x exp 



2tt m r 



N °9 (3.19) 
n(N,M) = c N H / dm r m; 3/2 e^ m ^ M -5(^J2 m r- M ">- 



r=1 M 
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Performing all the necessary steps, one recovers (3.12) with 

/(0 = 



0FErfi(* 1/4 ) - v^Erfi(e 1/4 ) - + ' " 



where Ern(z) = Erf (iz) j% is the imaginary error function, 



(3.20) 



/ 4tt 2 M 2 \ _ 4tt2 



and 



4vr 2 M 2 



(3.21) 



v^Mp 



(3.22) 



Therefore one has, 



S = Sf4V (3-23) 



3/ 2 V 4 



The constant 7 will depend on Mo, the lowest value of the mass for which the 
density of levels, p2{m r ) ) in (3.6) is valid, and on the constant, c. We have taken 
c = sJ M p , Mo = M p and found 7 ~ 0.089. 7 is found to decrease sharply with 
increasing Mq/M p . This behavior contrasts with the expected value of 717/3 ~ 1. 
The reason for this discrepancy may be traced to our use of the asymptotic level 
density. Recall that m 2 . < M 2 /N and we found, quite generally, that the density 
of states was maximized when N ~ M 2 /M 2 . This implies that m r ~ M p , which 
contradicts our use of (3.6). 

When only the lowest levels are occupied, we can approximate po{ m r) by a 
constant, say k/M p . This gives /(£) = C, 1 ^ 2 , b = k 2 and a = 1, and inserting these 
values into (3.17) and (3.18) yields 

s = s^(t)- < 3 ' 24 > 



To recover Hawking's temperature we must take k ~ y/87re. 
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4. Discussion 

We have shown that the canonical reduction of spherical geometries due to 
Kuchar 19 , combined with the introduction of an extrinsic time variable via the 
coupling to non- rotating dust as proposed by Kuchar and Brown 20 , leads to a re- 
markably simple description of the eternal Schwarzschild black hole. The Wheeler- 
DeWitt wave functional was seen to be expressible as a direct product state, which 
was interpreted as a collection of iV oscillator hierarchies, each of which could be 
described by a free, massless, complex scalar propagating in a "flat" two dimen- 
sional background (the internal, "metric" space) and confined to the interior of the 
Kruskal manifold. The canonical quantization program does not allow for an esti- 
mate of N, but, using techniques from the study of the statistical behavior of dual 
models, the Veneziano model and the statistical bootstrap model by Frautschi 16 
and Carlitz 17 , we discovered that the value of N that maximizes the density of 
states is given by N ~ M 2 /M 2 , i.e., the number of (lattice) sites is proportional to 
the area of the horizon, implying that all the degrees of freedom may be thought 
of as residing on the horizon itself. It is a remarkable result that probably has a 
deeper meaning. One is tempted to speculate, for example, that this, or some sim- 
ilar mechanism, may be a consequence of, or indeed a justification for, some form 
of a "holographic" principle. Moreover, because the total energy, M, of the black 
hole, is divided between these N lattice sites according to the dispersion relation 
X] r m r = M 2 , each oscillator is virtually in its ground state, with a small associ- 
ated degeneracy. The picture that emerges thus coincides also with Bekenstein's 
original way 3 of estimating the entropy by dividing the horizon into cells of Planck 
area, each of which has a small number of associated states. 

We were able to calculate the statistical entropy of the black hole by evaluating 
the microcanonical density of states of the system. Thus we recovered the area law 
of black hole thermodynamics, which is seen to be the consequence of the dispersion 
relation in (3.5) and the area quantization rule. It is noteworthy that all the states 
are explicitly displayed, both the mass quantization rules and the statistical entropy 
are recovered, and no explicit appeal to boundary states has had to be made in 
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this approach. 
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